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ABSTRACT 

We report on a sub-stellar companion search utilizing interferometric fringe-tracking 
astrometry acquired with Fine Guidance Sensor 3 (FGS 3) on the Hubble Space Telescope. 
Our targets were Proxima Centauri and Barnard's Star. We obtain absolute parallax 
values for Proxima Cen Tiabs = 0'.'7687 ± 0'.'0003 and for Barnard's Star T^ats = 0'.'5454 ± 
0'.'0003. 

Once low-amplitude instrumental systematic errors are identified and removed, our 
companion detection sensitivity is less than or equal to one Jupiter mass for periods 
longer than 60 days for Proxima Cen. Between the astrometry and the [Kiirster et al.| 



1999| radial velocity results we exclude all companions with M > O.SMjup for the range 



of periods 1 < P < 1000 days. For Barnard's Star our companion detection sensitivity 
is less than or equal to one Jupiter mass for periods longer than 150 days. Our null 
results for Barnard's Star are consistent with those of Gatewood (1995). 

Subject headings: astrometry — stars: individual (Proxima Centauri, Barnard's Star) 
— stars: late-type — stars: distances — stars: sub-stellar companions 
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1. Introduction 

Currently accepted theories predict that planetary 
syste ms are a natural b y-product of the formation of 
stars ( Lissauer, 1993 ). 



Black, 1995 reviews the im- 



portance of searches for extrasolar planets to theo- 
ries of solar system formation and discusses the lack 
of results to mid -1995. These search es h ave recently 
succeeded (e.g. Marcy et al., 1997 and Cochran et 



il., 1997). Radial velocity methods have been used 



for these detections. However, the derived planetary 
(or brown dwarf) masses are lower limits because of 
the unknown orbital inclination. These detections in- 
clude planets with minimum masses ranging from half 
a Jupiter mass (51 Peg, Mayor fc Queloz, 1995) to 



more than seven times the mass of Jupiter (70 Vir, 
Marry Rutler, 1 flSE) All nf these nhierts nrbit 



stars nf snlar spectral type^ pilthnugh they meiy tend 



to unusual metal richness (Gonzalez, 1996, Gonzalez 



1998 



1998). Recently Delfosse et al, 1998 and Marcy et al 



detected a planetary mass companion to the M4 
dwarf star Gl 876. Direct imaging has yielded two rel- 
atively high-mass objects, the 20 — 50 Afjupitcr brown 
dwarf, Gl 229B (Nakajima et al. 1995) and from Re- 



bolo et al. (1998) G 196-3B (15 
orb it M stars . 



40 M 



Jupiter J 



Both 



rhxt 



he relative position information provided by as- 
trometry resolves the mass uncertainty associated 
with radial velocity detections. An astrometric search 
recently evidenced a long-period Jupiter-mass com- 



panion to an M2 star, Lalande 21185 ( Gatewood 
|1996 ). Past examples of the successful application 



of astrometry to the discovery of stellar-mass unseen 
companions include Harrington (1977) and Lippincott 
(1977). 

Black and Scargle (1982) were the first to point 
out that Jupiter-like planets orbiting M stars might 
have short periods. They argued that scaling down 
the pre-planetary nebula thought to accompany the 
formation of stars would result in a gas-giant forming 
relatively close to a low-mass star. Such an object 
would have a period far shorter than that of Jupiter 
around the Sun. This was the motivation for our 
Hubble Space Telescope (HST) planet search program. 
The intervening years have changed this picture to 
include the formation of gas-giant planet cores well 



away from the parent star, even for M stars (Boss 



L995), with subsequent inward orbital migration (Lin 
al., 1996D. 



We present astrometry of Proxima Cen and Barnard's 



Star, including results of our astrometric searches for 
brown dwarf and planetary-mass companions. Those 
interested can trace the history of data acquisition 
and analysis over the past six years through a series of 
progress reports (Benedict et al, 1994a, 1994b, 1995, 
1997, and 1998c). 

Given the null result of an astrometric companion 
search carried out by Gatewood (1995) for Barnard's 
Star, the thirty four observation sets of this field 
nonetheless have significant value. Any systematics 
introduced by HST and/or FGS 3 should be present 
in the data for both Barnard's Star and Proxima Cen. 
In this sense Barnard's Star is a control in the exper- 
iment to detect low-mass companions orbiting Prox- 
ima Cen. 



Past direct detection observations include Van Bu 



ren et al. 1998, who observed Barnard's Star at lOyU 
and established companion upper limits 70 — SOMjup 
for separations 4 < a < 18 AU. A possible companion 
to Proxima Cen dete cted by the HST FG S used in a 
coronagraphic mode ( Schultz et al. 19*9^ ) is not con- 
firmed by our astrometric data, as was pointed out in 
that paper. Such a companion must possess an excep- 
tionally low mass to luminosity ratio. A subsequent 
study using HST WFPC-2 ( [Golimowski fc Schroeder 



1998| ) also failed to confirm this companion or find 



any as bright with separations 0.26 < a < 1.11 AU 
(0'.'2 <a< 0'.'85). 

In §2 we briefly review the astrometer, FGS 3, and 
discuss the data sets. We outline our calibration pro- 
cedures (§3), and present in §4 the astrometric model- 
ing of the reference frames. In §5 we derive parallaxes 
and proper motions for Proxima Cen and Barnard's 
Star. In §6 we explore the astrometric residuals for 
periodic phenomena indicative of companions. In §7 
we discuss the astrometric results, first the differences 
between HST and HIPPARCOS, and then our mass 
detection limits that rule out at the 0.1% level Brown 
Dwarf companions (M < 0.013Msun) to these two 
stars. See Opcnheimcr et al., 1999 for a discussion of 



the mass limits defining Brown Dwarfs. 

Tables ^ and ^ provide aliases and physical pa- 
rameters for our two science targets. The masses are 
based on the parallax results of this paper (§5). From 
those distances we obtain the V-band absolute mag- 
nitudes. My in Tables |l| and |^, then obtain masses 
from the recent lower main seque nce mass-luminosity 
relationsh ip of Henry et al. 1999 . Radii are from the 
models of Burrows et al. 1993, confirme d by the CM 
Draconis results of Metcalfe et al. 1996, 
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We time-tag our data with a modified Julian Date, 



an astrometric reference frame (see Benedict et al. 



JD ^ JD - 2444000, and abbreviate milli-second 



1998a for a more detailed discussion), we centroid 



of arc, mas, throughout. 



The Astrometer and the Data 



Our goal, 1 mas precision small-field astrometry. 



has^ been achieved, but not without Hlgiiiflcaut chal- 
lenges. Our observations were obtained with Fine 
Guidance Sensor 3 (FGS 3), a two-axis, white-light 
interferometer aboard HST. Bradley et al., 1991 pro- 



vide an overview of the FGS 3 instrument and Bene- 



dict et al., 1994a describe the astrometric capabilities 



of FGS 3 and typical data acquisition strategies. 

The coverage for both targets suffers from extended 
gaps, due to HST pointing constraints (described in 
Benedict et al., 1993) and other scheduling difficul- 
ties. For Proxima Gen the data now include 152 
shorter exposures secured over 4 years (March 1992 
to October 1997) and 15 longer exposures (July 1995 
to July 1996). Each orbit contains from two to four 
exposures. The longest exposure times pertain only 
to Proxima Gen observations obtained within Gon- 
tinuous Viewing Zone (CVZ) orbits. These specially 
scheduled orbits permit ^ 90 minutes on field, during 
which Proxima Gen was not occulted by the Earth. 
See Benedict et al. 1998 (Appendix 1, Table Al) for 
times of observation and exposure times for all Prox- 
ima Gen astrometry. A total of fifty-nine orbits have 
astrometric value. However, data acquired prior to 
mJD 8988 is of overall lower quality. During this era 
we had no on-orbit instrumetal astrometric calibra- 
tion and no independent scale assessment (§3.1), no 
assessment of the time- varying component of our as- 
trometric calibrations (§3.2), and employed only rudi- 
mentary (constant rate) intra-orbit drift correction 
(§3.3). 

Barnard's Star was monitored for three years (Febru- 
ary 1993 to April 1996), and observed three times 
during each of 35 orbits. Thirty-four of these orbits 
have astrometric value. Exposures range between 24 
and 123 seconds duration. See Benedict et al. 1998 
(Appendix 1, Table A2) for times of observation and 
exposure times for all Barnard's Star astrometry. 

3. Data Reduction and CaUbration Proce- 
dures 

Our data reduction and calibration procedures 
have evolved since the preliminary description given 



raw data, removing intra-observation jitter; apply an 
Optical Field Angle Distortion (OFAD) calibration 
( [Whipple et al., 1995| , [McArthur et al., 1997|); apply 



time dependent corrections to the OFAD (Mc Arthur 



in Benedict et al., 1994a. To remove systematics from 



et al., 1997); and correct for drift during each obser- 
vation set (intra-orbit drift). Additionally, we apply 
a lateral color correction (depending on star color in- 
dex) during the orbit-to-orbit astrometric modeling. 

3.1. The OFAD 

Until recently, no calibration star field with cat- 
aloged 1 mas precision astrometry, our desired per- 
formance goal, existed. We used FGS 3 to cafibrate 
itself with multiple observations of a distant star field 
(M35). A distant field was required so that dur- 
ing the 2-day duration of data acquisition, star po- 
sitions would not change. We obtained these data in 
early 1993 and reduced them with overlapping plate 
techniques to solve for distortion coefficients and star 
positions simultaneously. An astrometric catalog of 
27 stars with 1 mas positions, covering an area of 
1600 X 700 arcsec now exists. 

The aberration of the optical telescope assembly, 
along with the optics of the FGS, comprise the OFAD. 
The largest component of the design distortion, which 
consists of several arcseconds, is an effect that mimics 
a change in plate scale. The magnitude of non-linear, 
low frequency distortions is on the order of 0.5 sec- 
onds of arc over the FGS field of view. The OFAD 
is the most significant source of systematic error in 
fringe-tracking astrometry done with the FGS. We 
have adopted a pre-launch functional form originally 
developed by Perkin-Elmer, the builders of the FGS. 
It can be described (and modeled to the level of one 
millisecond of arc) by the two dimensional fifth order 
polynomial: 

X' = aoo + aioA + aoiF + aaoA^ + 002^^^ + anXY (1) 
+a3oX{X^ + Y^) + a2iX{X^ - Y^) 
+ai2YiY^ - X^) + aazY{Y^ + X^) 
-fa5oA(A2 + Y^f + aiiY{Y^ + X^ 
+a32A(A4 - Y^) + a2:iY{Y^ - X^) 
+auX{X^ - Y^ f + ao5Y{Y^ - X^ f 

Y' = boo + bioX + boiY + 620^' + bo2Y^ + bnXY (2) 
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+bi2Y{{Y^ - X^) + bo3Y{Y^ + X^) 

+ &50X(X2 ^ y2)2 ^ b4lY{Y^ + X^f 
+ b32X{iX^ - Y^) + b23Y(Y^ - X^) 

+buX{X^ - Y^ f + bo5Y{Y^ - X^ f 



where X, Y are the observed position within the FGS 
field of regard, X' ,Y' are the corrected position, and 
the numerical values of the coefhcients and bij are 
determined by calibration. The calibration observa- 
tions required 19 orbits and produced over 570 star 
position observations. We then employed GaussFit 



(Jefferys, 1988), a least squares and robust estima- 
tion package, to simultaneously estimate the terms 
in the OFAD equations and the star positions within 
M35. 

3.2. Maintaining the OFAD 

The FGS 3 graphite-epoxy optical bench was ex- 
pected to outgas for a period of time after the launch 
of HST. This outgassing was predicted to change the 
relative positions of optical components on the opti- 
cal bench. The result of whatever changes were tak- 
ing place was a change in scale. The amount of scale 
change was far too large to be due to true magnifica- 
tion changes in the HST optical assembly. To track 
these scale- like changes we revisit the M35 calibration 
field periodically, the ongoing LTSTAB (Long-Term 
STABility) series. LTSTABs will be required as long 
as it is desirable to do 1 mas precision astrometry with 
FGS 3. The result of this series is to model and re- 
move the slowly varying component of the OFAD, so 
that uncorrected distortions remain below 2 mas for 
the center of the FGS 3 field of regard. The LTSTAB 
series is also the diagnostic for deciding whether or 
not a new OFAD is required. 

3.3. Intra-orbit Drift Corrections 

A major improvement implemented since the ear- 
lier report ( Benedict ct al., 1994a ) is application of 
inter-observation (intra-orbit) drift corrections. Over 
the span of an orbit the positions reported by FGS 
1 and FGS 2 for the guide stars change. Benedict 



3t al., 1998a shows an observation set with drift ex- 



reducing the time available within an orbit to observe 
the science target. An observation set must contain 
visits to one or more astrometric reference stars, mul- 
tiple times during each observation sequence. Pre- 
suming no motion intrinsic to these stars over a span 
of 40 minutes, one determines drift and corrects the 
reference frame and target star for a drift generally 
quadratic with time. As a result we minimize the 
error budget contribution from drift to typically less 
than 1 mas. 

3.4. Lateral Color Calibration 

Since this calibration had not previously been ap- 
plied to our data, we describe it in some detail. Be- 
cause each FGS contains refractive elements (Bradley 
et al. 1991), it is possible that the position measured 
for a star could depend on its intrinsic color. Changes 
in position would depend on star color, but the direc- 
tion of shift is expected to be constant, relative to 
the FGS axes. This lateral color shift would be unim- 
portant, as long as target and reference stars have 
similar color. However, this is certainly not the case 
for the very red stars, Proxima Gen and Barnard's 
Star, hence our interest. Pre-launch ground testing 



ceeding 30 mas over the span of 36 minutes. The X 
and Y drift rates are generally dissimilar and usually 
not constant. The solution to this problem is sim- 
ple and effective, but imposes additional overhead. 



( Abramowicz-Reed, 1994) indicated for FGS 3 a lat- 
eral color effect predominantly in the X direction, 
with magnitude approximately 1 mas per unit change 
in B-V color index. 

An on-orbit test was designed and conducted in 
December 1991. Due to excessive spacecraft jitter 
(from a combination of the original solar arrays and 
insufficient damping in the telescope pointing control 
system) and insufficient knowledge of the OFAD, the 
results were inconclusive. We repeated the test in 
December 1994. Analyses of these data have yielded 
a lateral color calibration. 



3.4.1. The Lateral Color Field and Data 

Photometric and positional data are presented in 
Table ^ for the four stars in our chosen lateral color 
test field. We obtained POS mode data during three 
consecutive orbits on mJD 9708 (22 Dec 1994), plac- 
ing the four star asterism near each end and in the 
middle of FGS 3. We observed these additional posi- 
tions to explore Lateral Golor variations within FGS 
3. At each position HST was rotated -h60° and -120° 
from nominal roll. These large rolls were permitted, 
since on this date the field (on the ecliptic) was at 
the anti-sun. Of the nine data sets, one was so badly 
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affected by spacecraft motion to be unusable. 

3.4-2. The Lateral Color Model and Results 

Unfortunately, most of the originally proposed ref- 
erence stars proved too faint and were not observed. 
The test resulted in enough observations for a single 
overlapping plate reduction. The eight usable data 
sets were modeled by 

X' = X + ctx * {B ~ V) (3) 

Y' = Y + ctytc {B - V) (4) 

C = cos{R) x= X' - sin{R) * F' + c (5) 

ry = sin{R) * X' + cos{R) *Y + f (6) 

solving for R, spacecraft roll; c and f, offsets in X 
and Y; and global lateral color terms, ctx and cty, 
using the overlapping plate techniques described in 



Benedict et al., 1994a. Again using GaussFit we find 
ctx = -0.9 ± 0.2 and cty = -0.2 ± 0.3 mas. 
Figure |l| presents histograms of the residuals ob- 
tained by applying the model (equations 1 - 4) to the 
lateral color data. To demonstrate that there were 
no large remaining unmodeled systematic effects in 
these data we inspected the X and Y residuals plot- 
ted against X and Y position within FGS 3. 

4. The Barnard's Star and Proxima Cen As- 
trometric Reference Frames 

Figure || shows the distribution in FGS 3 pickle 
coordinates of the thirty-four sets of reference star 
measurements for the Barnard's Star reference frame. 
Figure ^ contains the corresponding distribution for 
the Proxima Cen reference frame. The circular pat- 
terns are impressed by the requirement that HST roll 
to keep its solar panels fully illuminated throughout 
the year. 

4.1. The Model 

From these data we determine the scale and rota- 
tion relative to the constraint plate for each observa- 
tion set within a single orbit. Since for Proxima Cen 
these observation sets span over five years (and over 
three years for Barnard's Star), we also include the 
effects of reference star parallax and proper motion. 

Our present model 



■q = d* X + e*Y + f - Py *TT - ^ly *t (8) 

differs slightly from that used previously on the 
Proxima Cen reference frame (Benedict et al., 1994a). 
X and Y are now corrected for lateral color (equations 
3 and 4) . The orientation of each data set is obtained 
from ground-based astrometry. Uncertainties in the 
field orientations are 0?03 and 0?09 for Proxima Cen 
and Barnard's Star. We obtain the parallax factors, 
Px and Py from a JPL Earth orbit predictor (DE200, 
Blandish, 199C). Finally we constrain S/i = and 



Stt = for the entire reference frame. 

Below we discuss the results of this modeling for 
the Proxima Cen and Barnard's Star reference frames 
in parallel. 

4.2. Scale Stability 

As in Benedict (1994) we form a scale-like param- 
eter, S, where 



5*= (ae-M)i/2 



(9) 



(, ^ a* X + b*Y 



* t 



(7) 



and a,h,d,e are the coefficients determined in equa- 
tion 5 and 6. S for both the Barnard's Star and Prox- 
ima Cen reference frames is plotted against time in 
days in Figure IJ. Error bars are the internal error, 
(75, derived from aa.b.d,e- Since both sets have an ar- 
bitrary scale zero point, 5*0, a shift has been applied 
to bring agreement in the mean scale in the range 
9000 < mJD < 10300. 

The most obvious feature is coincident with our 
early 1993 OFAD determination, at which time we 
moved from a predicted (from ground-based measure- 
ments at the manufacturing facility) to a measured 
OFAD. Servicing missions (SMI and SM2) have the 
potential to disrupt the position of FGS 3 within HST 
and introduce new volatiles. Neither SMI nor SM2 
seems to have an identifiable effect in Figure ^. 

From the scatter proximate to each observation 
date, it is obvious that the scale is indeterminate 
at the 5 x 10~^ level. Each observation set is ob- 
tained with the primary science targets (Proxima Cen 
and Barnard's Star) positioned within 5 arcsec of the 
pickle center. Thus, the effect of this scale indeter- 
minancy on their calibrated positions is on order 0.3 
mas. 

Since a search for periodic astrometric phenom- 
ena is a primary goal of this investigation, the quasi- 
periodic behavior of S should be characterized with 
the same tools we will use on Proxima Cen and 
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Barnard's Star. Our periodogram tools are described 
in §6 and Appendix Al.l. Figure |^ contains a peri- 
odogram characterizing the temporal behavior of S. 
Note the large amount of low-frequency power for 
800 < / < 2000(i~^. If some uncorrected fraction 
of this small variation is detectable in the measure- 
ments of our primary science targets, we would ex- 
pect to see it only in Proxima Cen. Observations of 
Barnard's Star are all post- mJD 8988, the location 
of the major break in the scale trend, pre- and post- 
OFAD. 

4.3. Reference Frame Results 

Tables ^ and ^ provide identifications, magnitudes, 
colors and relative positions for all reference stars. 
The relative positions within each frame have been 
rotated to J2000 RA and Dec, at epoch 1996.42 for 
the Barnard's Star field, 1996.79 for the Proxima Cen 
field. The Proxima Cen reference frame star m agni- 
tudes and colors are from Benedict et al., 1994a, as is 



the position of the star relative to which all positions 
are referenced, star ID-15. The Barnard's Star refer- 
ence star magnitudes and colors were obtained using 
the 0.8m CCD Prime Focus camera ( |Claver 1992| ) 
at McDonald Observatory. Positions are referenced 
to star ID-30, whose position is from the Guide Star 
Selection System Catalog ( Lasker et al., 199"c| ). 

Having applied the model (equation 5 and 6) we 
form histograms of the residuals (Barnard's Star, Fig- 
ure H; Proxima Cen, Figure 0) . No observations have 
been removed as outliers, but the overlap was carried 
out using the robust estimation technique described 
in Jefferys (1988). From these histograms we con- 
clude that our per-observation precision is 0.8 mas in 
X and 1.1 mas in Y. 

4.4. Assessing Reference Frame Residuals 

The OFAD reduces as-built HST telescope and 
FGS 3 distortions with magnitude ^ 1" to below 2 
mas ( [Mc Arthur et al., 1997 ) over much of the FGS 
3 field of regard. From the histograms (Figures ^ 
and ^ above) we have obtained correction at the ~ 1 
mas level in the region available at all HST rolls (an 
inscribed circle centered on the pickle-shaped FGS 
field of regard). To determine if there might be un- 
modeled, but eventually correctable, systematic ef- 
fects at the 1 mas level, we plotted the Barnard's Star 
and Proxima Cen reference frame X and Y residu- 
als against a number of spacecraft, instrumental, and 



astronomical parameters. These included X, Y posi- 
tion within the pickle; radial distance from the pickle 
center; reference star V magnitude and B-V color; 
and epoch of observation. We saw no obvious trends, 
other than an expected increase in positional uncer- 
tainty with reference star magnitude. 

5. Modeling the Parallax and Proper Motion 
of Barnard's Star and Proxima Cen 

Once we have determined plate constants from ap- 
plying equations 7 and 8 to the reference frames, we 
apply the transformations to the Proxima Cen and 
Barnard's Star measurements and solve for relative 
parallax and proper motion. At this step the lat- 
eral color correction is differential. The correction is 
based on the color difference between the target and 
the average color of the reference frame. Our relative 
parallax and proper motion results are presented in 
Table 0. 

Every small-field astrometric technique requires 
the following step; a correction from relative to ab- 
solute parallax. This correction is required because 
the reference frame stars have an intrinsic parallax. 
Ideally, all reference stars would be more distant than 
the desired science target star parallax precision, a sit- 
uation rarely provided by nature. To obtain science 
target parallaxes precise at the 0.5 mas or better level 
requires this correction. Faint (hence, in general, dis- 
tant) reference st ars have smaller co rrections. CCD 
techniques (e.g., Monet et al. 199*^ ) provide an ex- 
tremely faint reference frame {V ~ 18), but a rela- 
tively small field of view. FGS 3 provides a large field 
of view and access to moderately faint reference stars 
{V > 16). 

We adopt the corrections discussed and presented 



in the Yal e Parallax Catalog ( [Van Altena, Lee fc Hof 
Seit 1995| , Section 3.2, Fig. 2, hereafter YPC95). En- 
tering YPC95, Fig. 2, with the Barnard's Star galac- 
tic latitude, I = 13?99 and average magnitude for the 
reference frame, < Vref >= 14.2, we obtain a correc- 
tion to absolute of O'.'OOIO. 

As a test of the validity of this correction we (Cor- 
nell and Chappell) have obtained B, V photometry 
of the Barnard's Star reference frame stars using the 
0.8m CCD Prime Focus camera ( |Claver 199'^ at Mc- 
Donald Observatory. Classification spectra were ob- 
tained by one of us (Lee) with the WIYN Obser- 
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vatoryp^ multiobject spectrograph (MOS/Hydra). 
These data were reduced and then classified (spec- 
tral type and luminosity class) on the MK system 
by Lee. We find that the colors, magnitudes, and 
spectral types yield a reference frame whose bulk dis- 
tance properties agree with the YPC95 prediction. 
Table ^ collects reference star identification num- 
bers; V magnitudes; spectral types and luminosity 
classes; i3 — T/color indices; absolute magnitudes. My, 
and intrinsic color indices, (B-V)o from the color 
vs spectral type and luminosity class listed in Lang 
(1992); and resulting color excesses, Eb-v- Adopt- 
ing R — 3.1 — Av/Eq_y we obtain the listed to- 
tal V-band absorptions, Ay and, finally, the distance 
moduli, m-M, distances in pc, and parallaxes listed 
in the last two columns. The average parallax of the 
reference frame is < tt > = 1.2 mas. The difference, 
0.2 mas, between this and the YPC95 value, 1.0 mas, 
provides an estimate of the error in the correction to 
absolute. 

For Proxima Cen (/ = -1?93 and < Vref >= 14.5) 
we obtain a correction, -1-1.0 mas . 

Applying these corrections results in the absolute 
parallaxes listed in Table 0. The final formal uncer- 
tainties include the estimated error in the correction 
to absolute, 0.2 mas, RSS-ed with the relative paral- 
lax errors. We also list absolute parallax and proper 



motion from the HIPPARCOS catalog (Perryman et 
al. 1997) and the weighted, average parallax from 



YPC95. 

5.1. The Secular Acceleration of Barnard's 
Star 

With a large proper motion and parallax (fi = 
10'.'368y~^, TT — 0'.'5454) and large negative radial ve- 
locity {RV = —106.8kms~^), Barnard's Star is ex- 
pected to evidence a secular acceleration in the dec- 
lination component of the proper motion. This per- 
spective effect is characterized by 



/i = -2mfj.TrRV X 10" 



(10) 



and yields for Barnard's Star /i = 0'.'0012j/r^^. Gate- 
wood (1995) detected and confirmed the magnitude 
of this effect with 6 years of data. To determine if our 
data are sensitive to this effect we modify equation 8 



^*The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and 
the National Optical Astronomy Observatories. 



thusly, 

rj = d * X + e *Y + f — Py * TT — fly * t — fly * t * t (11) 

and find fiy = 0'.'0012 ± 0'.'0004yr-^ a 3 - cr detection 
from 3.3 y of data. 

Proxima Cen, with ^ = 3'.'851yr-\Tr = 0'.'7687, 
and RV = -21A:tos"^ has a predicted /i = 0'.'00013yr-2, 
undetectable by HST. 

6. Astrometric Companion Detection 

Applying the proper motion and parallax model 
(equations 7 and 11 for Barnard's Star, equations 7 
and 8 for Proxima Cen) yields residuals for Barnard's 
Star and Proxima Cen. We wish to test for the ex- 
istence of periodic behavior in these residuals, as- 
trometric perturbations indicative of planetary-mass 
companions. 

6.1. Forming Normal Points and Associated 
Uncertainties 

Every Barnard's Star orbit contained at least three 
observations. Before subjecting the residuals to a 
period-finding procedure we average the multiple ob- 
servations obtained at each epoch (the science target 
observations collected during a single orbit) to form 
normal points. This step removes the high-frequency 
power for 1// < 0.014(i, the typical time interval be- 
tween intra-orbit observations. Each normal point 
has an associated uncertainty that is the standard 
deviation of the intra-orbit observations. Resulting 
normal points and uncertainties are shown in Figure 

I- 

For Proxima Cen only those orbits acquired for 
t > mJD 9924 contained multiple observations of the 
science target. However, prior to that epoch the fre- 
quency of observation was far greater. Consequently, 
we trade dense temporal coverage at earlier epochs 
for more realistic normal points and uncertainties by 
averaging over 3-10 days. This approach results in 
the 28 normal points plotted in Figure ^ 

To summarize, Barnard's Star and Proxima Cen 
(for t > mJD 9924) normal points are always formed 
from data acquired within one orbit. Proxima Cen 
normal points prior to mJD 9924 are formed from 
data acquired over multiple orbits. 

6.2. Searching for a Perturbation 

The most obvious perturbation search is visual in- 
spection of the residuals plotted against time. We 
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look for trends or obvious periodicities. The resid- 
uals for Barnard's Star (Figure ||) appear random. 
The residuals for Proxima Cen (Figure ||) do not. We 
do not have a seamless set of observations for Prox- 
ima Cen. Data prior to mJD 8988 had no OFAD, no 
OFAD maintenance (LTSTABS), and only rudimen- 
tary (constant rate) drift correction. If we remove 
those data prior to mJD 8988 and redo the entire 
analysis, the residual vs. time plot becomes far flat- 
ter. Additionally a weighted fit of a parabola to all 
the data in Figure H yields statistically insignificant 
curvature. 

We next search for more subtle variations. Since 
the HST observations of Proxima Centauri and Barnard's 
Star are unevenly distributed in time, traditional 
Fourier transform techniques which assume regularly 
spaced data are inappropriate for period searching. 
The usual tactic of clamping the data to zero during 
large gaps in the sampling often leads to excessive 
power on scales of the lengths of the gaps (Press et 
al. 1992). Also, if we interpolated the data onto a 
regularly spaced grid, information would be lost and 
the sought signal effectively degraded. 

Schuster (1905) introduced the periodogram as a 
simple way of handling data sets with non-uniform 
sampling such that each data point receives equal 
weight. This approach, described in Appendix A 1.1, 
was used in a successful search for photometric vari- 



ability of Proxima Cen and Barnard's Star (Benedict 
et 41- 1998bD . 

In Figures ^ and ^l|, normalized periodograms 
are presented for the positional data (the normal 
points in Figures || and ^ for Barnard's Star and 
Proxima Cen. Superimposed on each spectrum are 
horizontal lines of constant false-positive rates com- 
puted using Equation A-4. Barnard's Star evidences 
large power at a frequency equivalent to one-half year 
and Proxima Cen has a significant signature at the 
lowest frequencies examined for perturbations. 

Before further exploring the causes for the strong 
signatures in the periodograms, we wish to see if these 
signatures persist for an alternate spectral analysis 
technique. A Bayesian spectrum analysis technique 
developed by Bretthorst (1988) and Jaynes (1987) 
uses probability theory to approach the problem of 
signal detection (see Appendix A1.2). This approach 
is most beneficial where the signal-to-noise ratio is 
small. Figure ^ contains the resulting Bayesian 
spectral estimate for Barnard's Star. In this case the 
false-positive levels are estimated using the Monte- 



Carlo method discussed in Appendix A1.2. Noting 
that the Bayesian and normalized periodogram ap- 
proaches produce similar results, we will use normal- 
ized periodograms in the rest of this paper. 

6.3. Significant Peaks in the Power Spectra 

Most of the power in the Barnard's Star peri- 
odogram occurs at / ~ O.by^^. The Barnard's Star 
field is unique in one respect, motion across the pickle. 
It is our operational constraint to place the science 
target always within a few arcsec of the pickle cen- 
ter. Barnard's Star has a proper motion in excess 
of 10"yr~^. This combination dragged the reference 
frame over 30" across the pickle during our monitor- 
ing. Consequently, the reference star positions were 
shifted by 30"across the pickle, while the Sun con- 
straint rotated the field several times through 360°. 
Thus, the plate constants derived will likely have pe- 
riodic OFAD deficiencies buried in them, causing the 
strong signal in the periodogram. 

Both parallax and spacecraft roll correlate with an 
integral fraction of one year. Figure |l3| shows x and 
y normal points plotted against spacecraft roll (roll 
modl80°). There is a clear systematic effect in Y, 
linear with roll. In X the effect is quadratic with 
roll. If we substitute the x and y residuals (CorrNP) 
to these linear fits as our new normal points, the re- 
sulting power spectra (Figure |lj) no longer contain 
significant power at / ~ 0.5y~^. This empirical cor- 
rection renders our search insensitive to perturbations 
with 1/2 year periods (/ ^ 0.006day~^). 

With a proper motion roughly one-third that of 
Barnard's Star, the Proxima Cen periodogram might 
also be expected to have some signal at / ~ 0.0Q6day~^. 
However, it does not. The most significant peak 
in the Proxima Cen periodogram (Figure 11) oc- 
curs at the lowest frequency we are able to probe, 
/ ^ 2000(i^^. The periodogram is detecting the secu- 
lar drift in the residuals seen in Figure |^. Comparing 
the periodograms for scale (Figure |^) and the astrom- 
etry residuals (Figure pT|), the peak in the astrometric 
residuals periodogram has a frequency nearly identi- 
cal to the peak in the scale variation periodogram. 
The amplitude of the secular trend periodogram is 
much less, indicating possible incomplete correction 
of scale-like variations. 



8 



6.4. Detection Limits from Power Spectra 



7.2. Companion Mass Detection Limits 



Monte-Carlo simulations were used to determine 
angular perturbation limits for companions that could 
hide in these power spectra. For each signal pe- 
riod and amplitude, 2000 simulated data sets of a 
sinusoidal signal with additive uncorrelated Gaussian 
noise were generated. The noise variance and sam- 
pling pattern were chosen to match the reference star 
observations (§4.3). The probability that a spectral 
peak with amplitude Sa could have arisen from noise 
is given by false positive rate P/p(5a). We define 
a miss rate, PmissiSA) as the fraction of trials for 
which the largest spectral peak either falls below a 
given spectral amplitude, 5"^, or occurs at another 
frequency due to noise fluctuations. The derived per- 
turbation amplitudes are those in which 95% of the 
simulations produced power in excess of that corre- 
sponding to a 1% false-positive rate. The detection 
amplitudes represent a limit on the smallest under- 
lying signal which could give rise to a given spec- 
tral amplitude. We find that the detection ampli- 
tudes are effectively insensitive to perturbation fre- 
quency. For Barnard's Star we obtain 1.25 mas, and 
for Proxima Cen, 1.0 mas for perturbation frequencies 
/ < 1500d-i. 

7. Discussion of Astrometric Results 

7.1. The Parallax and Proper Motions of Prox- 
ima Cen and Barnard's Star 

As shown in Table ^ we find parallaxes for Prox- 
ima Cen and Barnard's Star that differ substantially 
from values found by HIPPARCOS. The YPC95 par- 
allax values are weighted averages of four and six 
independent determinations for Proxima Cen and 
Barnard's Star, respectively. All values and associ- 
ated errors are plotted in Figure showing good 
agreement between the HST and YPC95 parallaxes. 
The 3.6 and 3.9 mas differences between HST and 
HIPPARCOS parallaxes may be explained in part by 



HIPPARCOS zonal errors (e.g., Narayanan & Gould 
19991 ). 

The proper motion differences can, in part, be at- 
tributed to the HST values being relative to a local, 
not a global, reference frame. That these proper mo- 
tion differences are at the 1 mas level, demonstrates 
agreement between the HST and HIPPARCOS scales 
to within one part in 4000. Scale differences cannot 
explain the parallax differences. 



The only confirmed planetary mass companion 
found thus far for an M star has a s hort period 
(Gl 876, P ^ 61'', iDelfosse ct al., 1998| and |Marcy 
et al., 1998 ). This may be partly a selection ef- 



fect, since M stars have been studied extensively with 
high-precision radial velocity techniques for only a 
few years. Nonetheless it is interesting that Fig- 
ures |lj and |ll| show no significant power for fre- 
quencies 60 < / < 12days^^. Neither is there 
power at the stellar rotation rates {ProtBam ^ 130^*, 
Protprox = 83.5 ± 0.5*^) inferred from the FGS pho- 



tometry described in Benedict et al. 1998b| . This is an 



expected result, since the star spots, with 1-2% con- 
trast would not be expected to perturb the observed 
photocenters of Proxima Cen and Barnard's Star. 

Our detection limits are expressed in angular mea- 
surements. We require primary star masses and dis- 
tances to translate to mass detection limits. We esti- 
mate the masses of Proxima Cen and Barnard's Star 
from the mass-luminosity relationship of Henry et al. 
(1999). The masses are averages of a mass estimate 
from My derived from our parallax values (Table |^ 
and V, jLeggett 1992|) an d a B-V =f(My) relation- 
ship (Henry et al. 1999), where B-V is also from 



Leggett (1992). We adopt Mprox = O.IIMq and 



Mb am 

= O.16M0. 





We wish to estimate the mass 
of a planetary companion, Mp, from detectable val- 
ues of a' , the semi- major axis of perturbation orbit 
in arcsec. From 



(M* + Mp)P^ = a-^ 



aNh = aMp 
a'(") = a{AUy 



(12) 

(13) 
(14) 



with P in years, a and a in AU, and M in solar masses, 
we derive Mp able to produce the angular perturba- 
tions detectable at a 1% false-positive level (Barnard's 
Star, 1.25 mas; Proxima Cen, 1.0 mas). Our results 
for Barnard's Star and Proxima Cen are shown in 
Figures ^ and |l^. On the latter we show detection 
limits from the Proxima Cen radial velocity compan- 
ion search program of Kiirster et al. (1999), showing 
the complementarity of the two searches. 

While fairly extended in period coverage, our re- 
sults are less so when stated in orbit size, as shown on 
the top axes of Figures |l^ and |l^. How far from the 
primaries our techniques probe bears directly on com- 
parison with searches carried out by more direct tech- 
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niques, particularly camera observation of the neigh- 
borhoods of these stars. A Saturn-mass companion 
with P = 2000^^ would orbit Proxima Cen at a ~ 1.5 
AU (~ 1'.'2). A half- Jupiter mass companion orbiting 
Barnard's Star with P = lOOO'', would have a ~ 1.1 
AU (~ 0'.'58). There is satisfactory overlap for Prox- 



ima Cen with the camera observations of Golimowski 
& Schroeder 1998| , Lcincrt et al. 1997, and Schultz 
et i J. 1998|. For Barnard's St ar there remains a gap 



since 
AU. 



8. Conclusions 



Van Buren et al. 1998 reach only inward to 4 



1. FGS 3, a white light interferometer on HST 
used in a fringe-tracking (POS) mode, produces 1 mas 
precision astrometry. 

2. Attaining this precision requires multiple re- 
visits to reference stars during an orbit, necessary to 
remove short time-scale temporal positional drift. 

3. Longer time-scale instrumental changes are 
monitored and removed through periodic revisits to 
our primary astrometric calibration field, M35. 

4. These techniques yield parallaxes for Barnard's 
Star and Proxima Cen with precision better than 
0.4 mas. HST parallaxes differ from the HIPPAR- 
(705 determinations by 3.9 and 3.6 mas, respectively. 
The i/S'T determination agrees with the YPC95 value 
within the errors for Barnard's Star. The secular ac- 
celeration of Barnard's Star was determined at the 
3-tT level in 3.3 years of observation and agrees with 
past determinations. 

5. We have examined the Barnard's Star and Prox- 
ima Cen astrometric residuals for periodic perturba- 
tions due to planetary mass companions, using nor- 
malized periodograms and Bayesian spectral analysis. 

6. We find statistically significant power at one- 
half year for Barnard's Star. The Barnard's Star 
residuals correlate with spacecraft roll. We tenta- 
tively identify this systematic effect with reference 
frame motion of over 30" across the pickle during our 
monitoring. Our empirical correction for this sytem- 
atic error renders us insensitive to companions with 
periods P ~ O.S*'. The periodogram of the corrected 
data contains no significant power at any inspected 
frequencies. 

7. For Proxima Cen we find statistically significant 
power at / ~ 2000(i~^. The source of this power can 
be seen in the residuals plotted against time. Much 
of the non-linearity comes from lower quality data se- 



cured early in our program. We point out the similar- 
ity of this signature with that found in a periodogram 
of scale- like variations in FGS 3. Together they sup- 
port the probable identification of a systematic error, 
not a long-period, low-mass companion. 

8. After discarding the signals from systematic 
errors, Monte Carlo simulations indicate that a 1 
mas amplitude perturbation with a frequency range 
2000d"i > / > 30d"i would be found in the Prox- 
ima Cen residuals at a 1% false positive level, with a 
5% miss rate. Corresponding sensitivity for Barnard's 
Star is 1.25 mas. 

9. For Proxima Cen, assuming Mprox = O.IIMq, 
we detect no companions more massive than Jupiter 
with orbital periods SO'' < P < lOOO'^. For P > 400^^ 
our detection limit is less than the mass of Saturn. 
These null results confirm the direct observations of 



Leinert et al. 1997 and Golimowski & Schrocdcr 1998 



They weaken the interpretation of the FOS coron- 
agraphic observations of Bchultz et al. 1998. Be- 



tween the astrometry and the Kiirster ct al. 1999 



radial velocity results we exclude all companions with 
M > 0.8Mj„p for the range of periods 1 < P < 1000 
days. 

10. For Barnard's Star, assuming M Bam — 0.16Mq, 
we detect no companions more massive than twice 
Jupiter with orbital periods SO'' < P < ISO''. For 
P > ISO'' our detection limit is less than the mass of 
Jupiter. 
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APPENDIX- Perturbation Search Techniques 



Appendix A1.2 - A Bayesian Approach 



The standard periodogram when apphed to non- 
uniformly sampled data has some unfortunate statis- 
tical shortcomings, including a sampling-dependent 
probability distribution for pure noise samples and 
non-translational invariance. Both of these problems 
are cured by a modified version of the periodogram 
introduced by Scargle (1982). We used this modified 
periodogram to analyze the HST data and briefly re- 
view the method in the next section. Wc also used 
a technique called Bayesian spectrum analysis which 
is based on probability theory (see Bretthorst 1988). 
This method provides a way of estimating the peri- 
odogram defined by Equation 8, but is more robust 
for very noisy data sets. 

Appendix Al.l - The Normalized Periodogram 

Scargle (1982) introduced a new normalization for 
the periodogram which produces an exponential prob- 
ability distribution for random Gaussian noise inde- 
pendent of the sampling pattern. He also included 
a time delay, defined below, which solves the trans- 
lational invariance problem mentioned above. With 
these modifications, the "normalized" periodogram is: 



(A-1) 



1 f[E,'2/jCOsa;(ij-r)]2 [£i Vj smoj{tj - t)] 



+ 



2 1^ E,- cos2 witj - t) ' ^ . sin^ co{tj - r) 
where the delay, r, is defined by 

J2j sin 2(jjtj 



tan(2a;T) 



J2j cos 2ujtj 



With the inclusion of the time delay, the normalized 
periodogram becomes equivalent to least-squares fit- 
ting of the data to sine waves. 

The resulting exponential form of the probability 
distribution is convenient since the false alarm prob- 
ability takes a particularly simple form. For uncor- 
related Gaussian noise, the probability, Pfp, that no 
peaks in the power spectrum have amplitudes greater 
than, Sa, is 



Pfp = l-(l-e-*^) 



(A-3) 



where Ngs is the number of independent frequencies 
searched (Scargle 1982). 



A Bayesian spectrum analysis technique developed 
by Bretthorst (1988) and Jaynes (1987) starts with 
the hypothesis H that a periodic signal exists in the 
a data set D. The goal of signal detection is to find the 
probability P{H\D,r) that the hypothesis is correct 
given the data and any prior information / that might 
exist, e.g. noise statistics. Bayes' theorem relates this 
probability to other calculable probabilities as follows: 



P{H\D,I) 



P{H\I)P{D\H,I) 
P{D\I) 



(A-4) 



where P{H\I) is the prior probability of the hypothe- 
sis, e.g. if one hypothesis is more likely than another 
based on, say, theoretical grounds; P{D\I) is the prior 
probability of the data which for this analysis is just 
a normalization constant; and P{D\H, I) is the direct 
probability of the data (also called the likelihood func- 
tion) and is the probability that the data set could be 
produced given the hypothesis and the prior informa- 
tion. 

The hypothesis for the planet search is that a si- 
nusoidal signal defined by f{t) = Bi sin ojt + B2 cos cot 
exists in the data. The probability that this hypoth- 
esis is true may be written as P{oj, Bi, B2\D, I) = 
P{H\D,I). However, since only a detection is re- 
quired, we really only want the probability that a 
signal with frequency cu exists regardless of the am- 
plitude and phase, i.e. P{uj\D,T). BJ refer to the un- 
wanted parameters Bi and B2 as "nuisance" param- 
eters and demonstrate that they can be eliminated 
from the calculation by a simple integration: 

(A-2) P{u;\D,I)= J J P{io,Bi,B2\D,I)dBidB2 (A-5) 



where P(a;, Bi, B2I-D, /) = P{H\D,I) is given by 
Equation 12. They assume that Bi and B2 have flat 
probability distributions and include this assumption 
in the prior /. This method of reducing the dimen- 
sionality of the problem is of course very useful useful 
for models which contain a large number of parame- 
ters. 

Bretthorst and Jaynes show that since the differ- 
ence between the model and the data is ideally just 
the noise, the likelihood function depends on the sum 
of the squares of the difference between the data and 
the model function f{t). By substituting the model 
f{t) into the resulting likelihood function, expanding 
appropriate terms, and integrating over the nuisance 
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parameters, they show that the desired probabihty is 
related to the periodogram defined in Equation A-1 
by: 



P{lo\DJ) oc 



1 



2CV) 



(2-^)72 



(A-6) 



where y is the mean of the data. In the derivation of 
this expression it is assumed that only the shape of the 
noise probability distribution is known and not noise 
variance. One of the advantages of this technique is 
that it provides an internal estimate of the accuracy 
of its predictions; a feature which direct periodogram 
and Fourier transform methods lack. 

Since we are not aware of a corresponding analytic 
form for false positive rates for the Bayesian tech- 
nique, Monte-Carlo simulations were employed to ob- 
tain those plotted in Figures 



12, Numerical simula- 



tions were first used to look for frequency-dependent 
false positive variations of the Bayesian spectral esti- 
mate. This was done to insure that the non-uniform 
spacing of the data points could not interact with the 
noise to produce excessive power at any of the fre- 
quencies of interest. Groups of artificial data sets 
were created using the sampling times of the ob- 
servations and contained only uncorrelated Gaussian 
noise. The variance of the noise was chosen to be 
the estimated noise variance of each of the observa- 
tions. Bayesian spectra were computed for each arti- 
ficial data set over the range of frequencies of interest. 
Then, the frequency and amplitude of the largest peak 
was recorded for each spectra. In order to estimate 
the distribution of the largest peaks as a function of 
frequency, the frequencies of the peaks were collected 
into A^cff bins where N^s is the number of independent 
frequencies in the spectral region of interest. 



Figure [12| shows the results of this analysis for the 
sampling times of the Proxima Centauri positional ob- 
servations. 2000 artificial data sets and spectra were 
created and 40 bins were used. The bottom panel of 
the figure shows the spectral amplitude correspond- 
ing to false positive rates of 0.1%, 0.3%, and 1.0% 
over frequencies of 40~^ to 400~^days~^. In the top 
panel the relative number of peaks are shown over 
the same frequency range. Since the variations in 
this quantity are consistent with counting statistics, 
we conclude that the false positive rates are essen- 
tially frequency independent over range of frequen- 
cies of interest. Thus, the false positive rates for the 
Bayesian technique shown in Figures ^ were based 
on the probability distribution of peak amplitudes 



found over the entire spectral range of interest. 
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Table 1 

Proxima Cen = a Cen C = GJ 551 = V645 Gen HIP 70890 



Parameter 



Value 



Reference 



V 

B-V 

Mv 

Sp.T. 

-^^Prox 

-'^Prox 

-Rprox 



11.09±0.03 
1.90 ± 0.04 
15.60 ±0.1 

M5Ve 

O.IIMq ±0.01 

O.OO1L0 

O.14i?0 



Leggett 1992 



Leggett 1992 



from TT, this paper 



Gliese & Jahreiss. 1991 



Henry et al. 1999 



Liebcrt & Probst 1987 



Burrows et al. 1993 



Table 2 

Barnard's Star = GJ 699 = G140-24 = LHS 57 = HIP 87937 



Parameter 


Value 


V 


9.55±0.03 


B-V 


1.73 ± 0.04 


Mv 


13.3 ±0.1 


Sp.T. 


M4Ve 


MBarn 


O.16M0 ±0.01 




O.OO46L0 


-RBarn 


O.19i?0 



Reference 

Leggett 1992 
Leggett 1992 
from TT, this paper 



Kirkpatrick & McCarthy 1994 


Henry et al. 1999 






Henry & McCarthy, 1993 




Burrows et al. 1993 





Table 3 

Lateral Color Calibration Stars 



ID 


V 


B-V 


RA (°) 


Dec (°) 


1 


10.55 


1.92 


89.8712 


22.5804 


2 


10.42 


0.18 


89.8640 


22.5777 


6 


14.25 


0.81 


89.8861 


22.5775 


9 


14.03 


1.51 


89.8734 


22.5629 
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Table 4 

Barnard's Star Reference Frame 



Star ID V B-V C ctj ?? ct, 



2 ID-31 15.47 1.09 -146.5100 0.0009 34.1564 0.0012 

3 ID-32 14.81 0.59 -142.1551 0.0009 -3.6702 0.0013 

4 ID-33 14.64 1.01 -170.9684 0.0011 -52.0054 0.0016 
5* ID-30 13.17 1.05 0.0000 0.0008 0.0000 0.0009 
6 ID-34 14.09 1.35 -72.6553 0.0007 94.8155 0.0007 
8 ID-36 11.55 0.32 -37.0964 0.0006 90.0358 0.0007 



*RA, Dec = 269.4768708, +4.6805700 (J2000) 



Table 5 

Proxima Cen Reference Frame 



Star ID V B-V ^ (^^ r] a, 



2 ID-13 15.10 0.72 -18.4819 0.0008 -74.9461 0.0014 

3 ID-14 15.76 0.87 -3.8168 0.0017 -32.3059 0.0021 
4* ID-15 14.36 0.66 0.0000 0.0005 0.0000 0.0008 

5 ID-16 14.58 0.52 70.4847 0.0009 -113.6714 0.0017 

6 ID-17 15.32 0.73 15.6880 0.0018 -3.8986 0.0022 

7 ID-18 15.14 1.53 -3.6323 0.0007 27.5414 0.0011 
9 ID-19 14.30 1.05 131.3827 0.0008 63.0020 0.0013 



*RA, Dec = 217.435381, -62.691332 (J2000) 
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Table 6 

Barnard's Star Reference Frame Absolute Parallax 



Star 


ID 


V 


Sp.T 


B-V 


My 


(B-V)o 




Av 


m-M 


dist(pc) 


ttC) 


2 


ID-31 


15.47 


KOV 


1.09 


5.9 


0.81 


0.28 


0.88 


9.57 


547 


0.0018 


3 


ID-32 


14.81 


FOIII 


0.59 


1.2 


0.3 


0.29 


0.90 


13.61 


3480 


0.0003 


4 


ID-33 


14.64 


G9V 


1.01 


5.7 


0.77 


0.24 


0.74 


8.94 


437 


0.0023 


5 


ID-30 


13.17 


G9IV 


1.05 


3.1 


0.85 


0.20 


0.63 


10.07 


773 


0.0013 


6 


ID-34 


14.09 


G9III 


1.35 


0.75 


0.97 


0.38 


1.17 


13.34 


2714 


0.0004 


8 


ID-36 


11.55 


A9III 


0.32 


1.35 


0.26 


0.06 


0.19 


10.2 


1005 


0.0010 



17 



Table 7 
Parallax and Proper Motion 





Proxima Cen 


Barnard's Star 


HST study duration 


5.6 y 


3.3 y 


number of observation sets 


59 


34 


ref. stars <V> 


14.5 


14.2 


ref. stars < B — V > 


+0.87 


+0.85 


Relative Parallax 


767.7 ± 0.2 


544.2 ± 0.2 mas 


corr to absolute 


1.0 ± 0.2 


1.2 ± 0.2 mas 


HST Absolute Parallax 


768.7 ± 0.3 


545.4 ± 0.3 mas 


Yale Parallax Catalog (1995) 


769.8 ± 6.1 


545.6 ± 1.3 mas 


HIP Parallax 


772.3 ± 2.4 


549.3 ± 1.6 mas 


HST Proper Motion 


3851.7 ± 0.1 


10370.0 ± 0.3 mas y-^ 


in p. a. 


78?46 ± 0?03 


4?45 ± 0?1 


HST Secular Acceleration 




12 ± 4 mas y~^ 


HIP Proper Motion 


3852.9 ± 2.3 


10368.6 ± 2.1 mas y-^ 


in p. a. 


78?50 ±0?03 


4?42 ± 0?07 



18 




Fig. 1 . Histograms of Lateral Color calibration ob 
servation residuals to the model (equations 5 and 6) 
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Fig. 2. — Barnard's Star reference star observations in 
FGS 3 pickle coordinates. The symbol shape identifies 
each star (hsted in Table 4). 
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Fig. 3. — Proxima Cen reference star observations in 
FGS 3 pickle coordinates. The symbol shape identifies 
each star (hsted in Table 5). 
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JD- 2440000 



t t I 



OFAD SMI 



Fig. 4. — Comparison of scales derived from equation 

9 for Prox Con and Barnard's Star reference frames. 
Major events are labeled: OFAD, first full-up astro- 
metric calibration of FGS 3; SMI and SM2, servicing 
missions. 



Scale Periodogram 




frequency (day ) 



Fig. 5. — Power spectrum of scales derived from equa- 
tion 9 for Prox Cen and Barnard's Star reference 
frames. The horizontal lines are levels of false-positive 
probability (equation A-4). 
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Fig. 6. — Histograms of X and Y residuals obtained 

from modeling the Barnard's Star refcrenc;e frame 
with equations 7 and 11. Distributions are fit with 
gaussians characterized by a as indicated. 
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Fig. 7. — Histograms of X and Y residuals obtained 

from modeling the Proxima Cen referenee frame with 
equations 7 and 8. Distributions are fit with gaussians 
characterized by a as indicated. 
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Barnard's Star 
Normal Points 
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Fig. 8. — Residual normal points in X and Y for 
Barnard's Star plotted against time. 
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Fig. 9. — Residual normal points in X and Y 
Proxima Cen plotted against time. 
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Fig. 10. — Lomb-Scargle periodogram of the Figure g 
normal points for Barnard's Star. The sohd hne is 
the sum of the X and Y power spectra. Horizontal 
lines denote levels of false positive probability from 
equation A-4. There is a significant signal at / ~ 
182d-i. 
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frequency (day ') 



Fig. 11. — Lomb-Scargle periodogram of the Proxima 
Cen normal points in Figure ^ The solid line is the 
sum of the X and Y power spectra. Horizontal lines 
denote levels of false positive probability. Note the 
significant signal at / ~ ISOOd^^. 
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Fig. 12. — Bayesian periodogram of the Figure ^ 
normal points for Barnard's Star. The sohd line is 
the sum of the X and Y power spectra. Horizontal 
lines denote levels of false-positive probability from 
Monte-Carlo simulations described in §A1.2 of the 
Appendix. This approach also detects significant sig- 
nal at / - 182d~i. 
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Fig. 13. — Barnard's Star normal points plotted 
against HST roll, modulo 180°. There are obvious 
trends, linear in Y, quadratic in X. The residuals 
(CorrNP) are adopted as the corrected normal points. 
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Barnard's Star 
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Fig. 14. — Lomb-Scargle periodogram of Barnard's 
Star (CorrNP) normal points corrected for systematic 
error depending on HST roll. The solid line is the sum 
of the X and Y power spectra. Horizontal lines denote 
levels of false positive probability. The correction has 
effectively removed the power at / ~ 182d^^ seen in 
Figure O. 



32 



HST HIPPARCOS YPC95 



0.774 
0.772 
0.770 - 
0.768 - 
^ 0.766- 

Oi 0.764 - 

O 

c3 "-^^^ " 

^ 0.550- 

0.548 - 
0.546 - 
0.544- 

0.542 - 
0.540 - 



Proxima Cen 
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Fig. 15. — Absolute parallaxes for Proxima Cen and 
Barnard's Star determined by HST (left) and HIP- 
PARCOS (center), compared with the YPC95 (right). 
The horizontal dashed lines are the weighted averages 
of the three independent determinations. 
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Monte-Carlo simulations indicate that 
an 1.25 mas amplitude perturbation 
would be found in the data at a 
1% false positive level, with a 5% miss rate. 
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Fig. 16. — Barnard's Star companion detection mass 
limits. The top axis gives the size of the semima- 
jor axis of a 1 Mj„p secondary. Monte-Carlo simula- 
tions indicate that an 1.25 mas amphtudc perturba- 
tion would be found in the data at a 1% false positive 
level, with a 5% miss rate. 
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Monte-Carlo simulations indicate that 
an 1 .0 mas amplitude perturbation 
would be found at a 1 % false positive 
level, with a 5% miss rate. 
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Fig. 17. — Proxima Cen companion detection mass 
limits. The top axis gives the size of the semima- 
jor axis of a 1 Mjup secondary. Monte-Carlo simula- 
tions indicate that an 1.0 mas amplitude perturbation 
would be found at a 1% false positive level, with a 5% 
miss rate. The dashed line presents the limits deter- 
mined by the radial velocity program of 
al. IT999I 
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